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^ ■ Abstract 

I The Hirota bilinear formalism and soliton solutions for a generalized 

■ Volterra system is presented. Also, starting from the soliton solutions, we 

. obtain a class of nonsingular rational solutions using the " long wave limit" 

procedure of Ablowitz and Satsuma, and appropriate "gauge" transfor- 
mations. Their properties are also discussed and it is shown that these 
solutions interact elastically with no phase shift. 



o 
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1 Introduction 



X 

■ It is known that with the development of sohton theory significant progress has 

been made in finding special solutions of integrable nonlinear evolution equations 
which include soliton solutions, rational solutions, similarity solutions and so on, 
using I. S. T.[p, Hirota bilinear formalism 0], Backlund transformations 0] and 
Wronskian formalism f^, [§. 

The class of rational solutions was firstly investigated for KdV equation 
Avery simple way to find them was developed by Ablowitz and Satsuma 0, and 
consist in taking the "long wave" limit in the multisoliton solution. This method 
was used succesfuUy to find the rational solutions also for other completely inte- 
grable systems. Although, for continuous case, many results have been obtained 
with respect to finding them, a relatively less work for discrete systems has been 
done §, 0. 

Unlike the (2+l)-dimensional continuous systems where the dynamics and 
phenomenology of rational solutions ("lump" solutions) are well studied 
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1^, the dynamics of (l+l)-dimensional nonsingular rational solutions was very 
little investigated |]12|. In this paper we'll derive the bilinear formalism for a 
(l+l)-dimensional discrete generalized Volterra system, which represent a cou- 



pled system of generalized Volterra equations proposed by M.Wadati [T^, [0. 
From its soliton solutions we'll obtain a class of real and complex rational non- 
singular solutions using Ablowitz-Satsuma limiting procedure. Here the crucial 
step is the existence of many nontrivial arbitrary constants and "gauge" trans- 
formations. It is shown that these solutions collide elastically with no phase shift 
using a "head-on" collision of two 1-rational nonsingular solutions. Also their real 
and complex structure is preserved. The properties, generalizations and possible 
developments are also discussed. 

The paper is organized as follows. In the section 2 is constructed the bilinear 
form of the system by a reduction to a variant of Ablowitz-Ladik system with 
nonzero boundary condition. In section 3 the construction of rational solutions 
is presented. In section 4 properties and dynamics are discussed and, finally, 
conclusions are given in section 5. 



2 Bilinear formalism 

The generalized Volterra system under consideration is: 

Qn = {a + bQn + cQl){Rn+i - Rn-l) 

Rn = {a + bRn + CRI) {Qn+1 - Qn-l) (1) 

with Qn, Rn ^ as n goes to infinity and a, b, c are real constants. 
For waves propagating in one direction we may put Qn = Rn and we'll find a 
generalized Volterra equation proposed firstly by M. Wadati [|13|, fT^ . 
By the following transformations: 



Un = ^j^^{Qn + b/2c) 
:{Rn + b/2c) 



V4:ac - 62 

Aac — 6^, 

t - ^(^^) (2) 



and putting ao = 

the system (1) becomes: 



Vn = {1+ f^) {Un+1 - Un-l) (3) 
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with Un,Vn — > tto as n goes to infinity. In order to have ao real we consider that 
c, b are nonzero and 4ac — is positive. The system (3) represent a variant of 
Ablowitz-Ladik system [15|, with nonvanishing boundary conditions. 
We are constructing the bilinear form of this system by taking: 

id f fn\ 

Un = ao- -— log(— ) 

2 dt Qn 

where we assume that fn/gn and f'nl Qn ~^ const, as n goes to infinity. 

Plugging (4) into (3), integrating with respect to t and taking into account 
the nonzero boundary conditions we'll get: 

tan-i un = tan-i ao - log (4^^^) 

^ 9n+lJn-l 
4.-1 4- ~1 ^1 / fn+l9n-l\ 

tan w„ = tan ao - - log ( — ) 

2 dn+lln-l 

Using (4) again we'll find after some algebra: 

DtfWn = 2(1 + ao) Sinh Dnfn9n 

(coshD^ + iaosmhDn)fng'n = fn9n 

(cosh Dn + iao Sinh Dn)fn9n = fWn (5) 

which represent the Hirota form of the system (3). We can see that for ao = 
we can recover the self-dual network form proposed by Hirota 0. 
The 1-soliton solution are given by: 

/„ = 1 + exp(?7i + 0i) 
= 1 + exp{r]i + (f)[) 
g^ = l + exp(?7i + ipi) 

g'^ = 1 + exp{r]i + ^[) (6) 



where 



rjj = 2ej (1 + ao) sinh Pjt + PjU + r/° 
ao ej + cosh Pj i 
2 ^ sinh Pj 2 ^ 

ao 1 + e,- cosh Pj ie, 
^-P^^ = T^-- sinhP, +f^-- 
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where the index j labels the number of the soliton i.e. for 1-soliton solution j — 1, 
for 2-soliton j = 1,2 and so on and Sj = ±1 indicates the propagation direction. 
On the other hand we've choosen the phases such that any multisoliton solution 
to be real. Kj represents arbitrary real constants which does not depend on n or 
t. The procedure of finding rational solutions relies mainly on this arbitrariness. 
Now we can proceed to the 2-soliton solution which has the form: 



cxp(?7i + 0i) + exp{r]2 + 02) + exp(?7i + r]2 + (f^i + h + ^12) 
cxp(?7i + 0i) + exp{r]2 + (P2) + exp(?7i + 772 + 0'i + 02 + A12) 
exp(?7i + ^pi) + exp{r]2 + ^2) + exp(r7i + 772 + ^"1 + V^2 + ^12) 



fn=l 

gn = i 

= 1 + exp(?7i + il)'^) + exp(?72 + ^"2) + exp(?]i + 772 + + ^2 + ^ 



12 



where 



sinh 



sinh 



P1-P2 ' 

. 2 ' 
■ Pi +P2 ' 



exp Ax2 

for solitons propagating in the same directions (eie2 

"cosh(^^) 



1) and: 



exp A 
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cosh(^^) 



-1^ 



(8) 



for solitons propagating in opposite directions (ei €2 

Using the classical procedure we can construct any N-soliton solution but for 
our purposes 1-soliton and 2-soliton solutions are suficient. 



3 Rational Solutions 

The "long wave limit" method of Ablowitz and Satsuma we'll use henceforth 
consists in three main steps. In order to obtain a N-rational solution from N- 
soliton solution we have to: 

• expand the Hirota form of the N-soliton solution in power series of Pj with 
j = l,...,iV 

• choose the phase constants such that all 0{P^^ with k less than N to be 
zero. 

• making Pj ^ we recover the N-rational solution from the 0{Pj) terms 
with / > TV 

We arc focusing on /„ and Qn- The procedure for //^ and (7^ works in the same 
manner. Thus, for 1-rational solutions we'll expand /„ and gf„ given by (6). We'll 
use the following notations: 



9j = 2ej{l + al)t + n 
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= exp77° 
Pj = % 

and the expansions will be in power of 5 and pj are 0(1) terms. 
So, up to 0{5): 



fn^i + Ci{i + spiei + o{s^)) 



2pi 2 12 

where (*) means changing the sign of the terms multiplicated with i (complex 
conjugate as if Q would be real for j = 1, 2). 

For ei = 1 we must takeXi proportional with S. To cancel 0(1) terms we're 
taking Ci = — I/ckq. Thus when 5 — > we obtain: 

id 9i+i/2ao 4ao(l + ao) /n\ 

2c«i b*! — z/2q;o AaQOf + 1 

In the same manner f„ = m„ so, up to a constant and time scaling factor, 

4ao(l + ao) 



Qn — Rn^ 



which is a nonsingular real rational solution. 

Forei = — 1 the situation is dramatically changed. In this case the expansions 
becomes: 

= 1 + '-K^Ci + \k^CiPi{Oi - iao/2)S + 0{S^) 

9n (/n) 

We have no divergent terms in 6 so we can take Ki ~ 0(1). But in this case 
the 0(1) terms appear to be complex conjugate in and gn- Accordingly there 
is no Ci to cancel both of them. To cope with this situation, we can see that Un 
and Vn are invariant under the following "gauge" transformation: 

k^k^ (10) 

Qn Qn 

where A is an arbitrary complex function which does not depend on t. 
Now, taking C,i — 2i/Ki we'll cancel the 0(1) terms only in /„. So, 

/n = -Pl(^l-^«o/2)5 + 0(<52) 

gn = 2 + Pi{ei+ia^/2)5 + 0{5'') 
Choosing K—1/5 and then (5^0 we'll find: 

fn ^ - in ^CtOx 
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and: 



"0 - ^^^°g(^i ~ ^ao/2) 



Vn = K)* (11) 

which represents complex sohtary waves. The physical significance of this type 
of solutions is not clear. The form of Qn and Rn will be given up to constants: 

n 1 + ^0 , , gi(l + Q:g) 

In order to study "head-on" coUision of the solutions (9) and (11) we'll con- 
struct the 2-rational solution from 2-soliton solution with solitons having opposite 
directions i. e. /„ and (?„ are given by (8) and phase factors by (7) with 6162 = — 1. 
Expanding and gn up to 0{6^) we'll take Ki ~ 0{6) in order to cut divergent 
behaviour of type in cxp (pi 

fn = Mo + MiS + M2S^ + 0{d^) 

Qn = (/n)* 

where: 

Mo = 1 + Ciao + ^i^2C2 + |i^2ClC2 

i 1 i 

Ml = -piCi + OiQpiCiOi + -aoC2K2p2 + -^iiiQiJ<^i + 

+ ^(«0?^ - l)PlClC2 + ^«0i^2ClC2(Pl^l +P2^2) (12) 

4 Pi 2 

The expresion ofM2 is very complicated. Anyway after choosing the phase con- 
stants it will be strongly simplified. 

Thus we can cancel 0(1) terms both in /„ and gn by taking C^i = —l/ao but 
we cannot do this for 0{d) terms. Taking K2 ~ 0(1) and (2 = 2i/-ft'2 we'll do it 
for fn only. Introducing a "gauge " transformation in the form A = 1/5 we will 
recover from M2 and gn from (Mi)*. So, the rational solution will be: 

/n = ^1^2-^ + ^(--«0^l) 

4 2 ctQ 



and 



9n = ei-— (13) 

2q;o 



f 

J n 



2q;o 
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4 Dynamics and Properties 



The solutions (13) and (14) represents a "head-on" colhsion of two rational soli- 
tary waves one being real and the other complex, having same but opposite 
velocities. To see this we're keeping 9i fixed in (13) and make t ±oo. We'll 
find that (13) becomes: 

i d ^ .9i + V2«0n 

which represent exactly the real solitary wave(^). 

Also, keeping 62 being fixed and making t — ^ ±00 we'll get: 

ao- log (^2 - ^ao/2) 

which means that we have no phase shift at this "head-on" collision so, we have 
a (l+l)-dimensional discrete "lump-type" dynamics. Also, the real and complex 
character is preserved. We don't know exactly if does exist real nonsingular 
rational solutions for both directions. 



5 Conclusions 

To summarize, we've showed that rational nonsingular solutions both real and 
complex exists even for discrete systems. So this system posses both solitons 
and nonsingular (l-(-l)-dimensional discrete rational solutions decaying to zero as 
0(n~^) and their character being real and complex, depending on the propagation 
direction. 

We've discussed the dynamics of this rational solutions and we've showed 
that we have no phase shift at the interaction, specific for the (2-|-l)-dimensional 
continuous "lumps" solutions. Unfortunately we don't know any condition of 
existence for this type of solutions and why they occur in (l-l-l)-dimensional sys- 
tems. Its continuum limit could be viewed as a coupled pair of mKdV equations 
with nonzero boundary conditions. We can see that the form of 1-rational so- 
lution (the real one) has the same structure as those of mKdV given by Ono 
[0 and Ablowitz-Satsuma 0, so the continuum does not affect it drastically. 
Probably a setting up of this type of solutions in a I. S. T. framework will clarify 
much more about their form, dynamics, stability. 
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